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Abstract 

We apply light-cone quantization to a 1-1-1 dimensional supersymmetric field theory 
of large N matrices. We provide some preliminary numerical evidence that when the 
coupling constant is tuned to a critical value, this model describes a 2-|-l dimensional 
non-critical superstring. 


Matrix models have been remarkably successful |l| in providing us with an understanding 
of non-critical bosonic string theories in dimensions less than or equal to two. Their extension 
to dimensions greater than two, however, runs into a number of difficulties, some of them 
related to the appearance of tachyons in the bosonic string spectrum. Quantization in the 
inhnite momentum frame is a promising approach to these models, where strings emerge as 
bound states of more elementary “string bits” Bi- Studies of the c = 2 matrix model using 
this technique have indeed revealed the expected tachyon problem 1§. 

While the bosonic strings are unstable above two dimensions, one expects this instability 
to be cured by target space supersymmetry. Full space-time supersymmetry requires that 
the string theory be ten-dimensional (with a possibility of subsequent compactification).[] It 
is well known, however, that in c-f 1 dimensional non-critical string theories the Liouville co¬ 
ordinate enters on a different footing from other dimensions, so that there is no full Poincare 
invariance. Therefore, one cannot demand a full c -f- 1 dimensional supersymmetry. The 
most one could require is the c dimensional “space supersymmetry” which was considered 
in ref. [0. While this work utilized the NSR approach to the super-Liouville theory, it is 
interesting to ask whether an equivalent Green-Schwartz formulation exists |^, An even 
more ambitious goal is to construct matrix models which describe non-critical superstring 
theories. In this paper we make a small step in this direction by formulating and numeri¬ 
cally studying a two-dimensional supersymmetric matrix held theory which is expected to 
describe 2-1-1 dimensional non-critical superstrings. 

Our model is a higher-dimensional generalization of the Marinari-Parisi model (the super- 
symmetric quantum mechanics of a large-A hermitian matrix) |^. It is not entirely clear 
whether the continuum limit of the Marinari-Parisi model describes a theory with target 
space supersymmetry; in fact, the conventional understanding is that this model provides an 
alternative description of the bosonic c = 0 string . There are some signs, however, of 

remaining target space supersymmetry [^. Since one-dimensional supersymmetry is quite 
trivial, we hope that going one dimension higher will reveal a more interesting structure. 


successful construction of critical superstrings from string bit models has been found recently 
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The price we have to pay is that the model is no longer exactly solvable. Luckily, using 
the techniques of light-cone quantization, we can still set up a scheme for calculating the 
spectrum of the matrix model in the large-limit. Our approach is largely an extension of 
the program of ii, with the advantage that the supersymmetry cures the sickness of the 
model. 

The c = 2 matrix model studied in ^ is dehned as a two-dimensional held theory with 
the euclidean action 


^ = / (TxTi 




where is an iV x iV hermitian matrix held. The connection of this model with 

triangulated random surfaces follows, as usual, after identifying the Feynman graphs with 
the graphs dual to triangulations. At the leading order in N, we obtain a sum over the 
planar triangulated random surfaces embedded in two dimensions. If the tadpole graphs are 
discarded, as they should be because they do not correspond to good triangulations, then 
the entire perturbative expansion is hnite. This is similar to what we hnd in the matrix 
models for c < 1. Therefore, it is sensible to look for a singularity in the sum over the planar 
graphs as A approaches some critical value Ac. For the c = 1 model the spectrum of energies 
becomes continuous at A = Ac. This is because we are really dealing with a 1 -|-1 dimensional 
string theory whose center-of-mass mode satishes |T^ 


e^ = pI 


where Pfp is the Liouville momentum which possesses continuous spectrum. By analogy one 
might expect that in the c = 2 model the spectrum of = 2P+P_ becomes continuous at 
A = Ac. Light-cone quantization of the model revealed, however, that the mass-squared of 
the lowest state becomes negative for sufficiently large A and, in fact, tends to —oo as A —>• Ac 
0. Although it is still possible that a continuum of states is forming around = —cx), it 
was difficult to draw definite conclusions from numerical studies of the spectrum. Further 
evidence for a phase transition at some critical value of the coupling was given [|T3 based 


on an approximate identihcation of the c = 2 model with a one-dimensional spin chain. 
However, the fact remains that the corresponding continuum theory, even if it exists in some 
sense, is a very sick theory containing a tachyon. 

In this article, we attempt to cure this sickness by considering a 1 -|-1 dimensional large-A^ 
hermitian matrix held theory with (1,1) supersymmetry 


[P+,P_] = [P±,Q±]={Q+,Q_} = 0, 

Ql = P+, Ql=P.. 

From the above algebra it follows that the operators P+ and P_ are positive dehnite, since 
they are squares of hermitian supercharges Q+ and Q_. This implies that = 2P+P_ is 
also positive dehnite. As we tune the parameters, we expect the theory to become ehectively 
2-1-1 dimensional due to the appearance of the Liouville mode which makes the spectrum 
of continuous. Furthermore, we expect the mass-squared of the lightest state to vanish 
at this critical point, since the appearance of new massless states is typically a signature 
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of the phase transition. Indeed, the continnons spectrnm with a vanishing mass gap was 
observed in the c = 1 [T^ and the Marinari-Parisi models 0. In this article, we present 
some evidence that snch massless continnnm of states does appear in the spectrnm of 1 + 1 
dimensional snpersymmetric matrix held theory. 

The simplest matrix held theory with the above snpersymmetry is given by the action 




(fxd^eii 


D^D^ + W{^) 


where $ is a (1,1) matrix snperheld 

‘hp = + O'^ij + dOFij, 

and Tjj is a matrix whose elements are two-component spinors 








In light-cone qnantization only T- tnrns ont to be dynamical. D is the covariant derivative 


where P^ are the two dimensional Dirac matrices in the Majorana representation 


po = 




P^ = 


We consider the simplest snperpotential W (d*) given by 


W(d>) = 


A 


3^/N 


$3 


The Feynman graphs which arise from the pertnrbative expansion of the path integral gener¬ 
ate triangnlated random snrfaces embedded in the snperspace with coordinates (x°, x^, 6^, 9“^). 
We expect that, as A is tnned to some critical valne Ac, the size of a typical graph diverges 
so that the world sheet continnnm limit may be dehned. The dimension two terms in the 
world sheet action are hxed by the target-space snpersymmetry (they may also be inferred 
from the snperspace propagator). 




world sheet — 


(faVhh°‘^ zi {daX^ — iOV^daO^ {d^x^ — iOV^dpO^ iz2da9dfj9 , ( 1 ) 


where Zi and Z 2 are normalization constants, and the metric is regarded as a dynamical 
variable. This action resembles the Green-Schwartz action. However, it has no Wess-Znmino 
term or manifest K-symmetry. Nevertheless, if we hnd anything non-trivial from onr model, 
it is a kind of non-critical Green-Schwartz snperstring. 

In terms of component helds, the matrix model action is written 


5 = / d^xTr 


\{9,4>r- + - \vp4>) - 


( 2 ) 
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where ^(0) = hh'(0) = ficj) — ■^0^- As such, it could also be thought of as a matrix 


analogue of the model considered in |]^ with zero central charge. The above Lagrangian 
contains cubic and quartic interaction terms which depend linearly and quadratically on the 
coupling constant A. One could now proceed to enumerate the states and determine the 
matrix elements of P_ as was done in theories without supersymmetry |^, Instead, we 
prefer to determine the matrix elements of supercharges for the enumerated states, as was 
recently proposed in |]T3. This has a number of advantages. Firstly, the supercurrent contains 


only terms up to cubic order in helds, and depends only linearly on the coupling constant 
A. Secondly, the supercharge is a more fundamental dynamical object, being the square 
root of the Hamiltonian. We hnd it particularly convenient to consider the combination 
/ = \/2iQj^Q_. This operator is bosonic, commutes simultaneously with and P_, and its 
square is equal to M^. I contains more information than because it hxes the ambiguity 
of sign when one takes the square root. When we search for evidence for the continuous 
spectrum, this fact will be rather useful. 

The action is invariant under the standard supersymmetry transformation 

^ ^ ij 

The supercurrent can easily be shown to be 

1 Tr + iV{4>)n] . 

v2 

The light-cone components of the supercharges are 

Q-= j dx- : Tr [^2 (d.cj)) T. 


Q+ = J dx- : Tr [H(0)T_] : 


Expanding in normal modes^ 




{x ) = 


foo dk_ 


Vo ySPT 


(aij{k-)e 


—ik—X I t 


+ a]i{k-)e 


ik—X 


1 T / 1 dk_ f 

) = -7^ L —m 




a / 2 ^ Vo a /2 


—ik—X 


+ bUk_)e 


ik—X 


the supercharges become 

iQ- = J dk'/k (blj{k)aij{k) — alj{k)bij{k)^ 


y 


^/N ■ 


dkidk2 


^The factor of 1/a/2 appears due to the unusual normalization of the kinetic term for in (|^). 
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y/kik2 


sjki{ki + k2) 

1 

k2{ki + ^ 2 ) 


a\j{ki)a]^{k2)hik{ki + ^ 2 ) + h\f.{ki + k2)aij{ki)ajk{k2 


a\j{ki)h]^{k2)aik{ki + /C 2 ) + + k2)aij{ki)hjk{k2 


hUki) a]^{k2)aik{ki + ^ 2 ) + + ^2)%(fci)aifc(fc2) 


We have dropped the subscript “minus” on k for brevity, set p = 1, and introduced the 
dimensionless coupling parameter y = A/(2y^p). The theory may be regulated in the 
infrared by compactifying the x~ direction and imposing periodic boundary conditions on 
the helds. The momentum k then takes on discrete values n/L, and the integrals above are 
replaced by sums. For massive particles we may discard the u = 0 modes, which give inhnite 
light-cone energy, and restrict to positive integer n. This regularization was also shown to 
preserve supersymmetry |^. We introduce discretized oscillators 


Ai^n) = -^=aij{k = n/L), Bij{n) = -^bij{k = n/L), 




which satisfy commutation relations 


[^Ip'(u), (u )] bnn'^ik^jl 

, B )j- bjin'bikbjl- 


(3) 


The discretized expressions for the supercharges are 
iQ- = XI {Bjj{n)Aij{n) - Alj{n)Bij{n)) 


n=l 


= X ^ (Bij{n)Aij{n) + Alj{n)Bij{n)) 

_ 1 \/Ti ^ ' 


n=l 


-i/X X 

m=i n2=l 

1 


y/nin2 L 
1 

^Jniini -y 72 , 2 ) 
1 

n 2 (ni -y 72,2) 


A\^{ni)A]^{n2)Bik{ni + U 2 ) -F B\^{ni + n2)Aij{ni)Ajk{n2) 


A\^{ni)B]f.{n2)Aik{ni + 712 ) -1- ^^(ni -F n2)Aij{ni)Bjk{n2) 
B\j{ni)A]f.{n2)Aik{ni + 712 ) -h ^^(ni -F n2)Bij{ni)Ajk{n2) 


The Hilbert space of this model can be constructed by acting on the Fock vacuum with 
mode creation operators A\j and B\y As in ref. @, ||, we restrict our attention to singlet 
states, [| i.e. states of the form 


1$) ~ Tr ( 711 )^^( 772 ) 5 '^(ns) • ■ ■ A^(7ii) |0). 


^The non-singlet states do not have an obvious stringy interpretation. As discussed in ref. Q they may 
be eliminated by adding some amount of the gauge interaction. 
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The light-cone momentum of such a state is given hy P- = K/L, where K = and 

we restrict ourselves to states of fixed P_. The integer K plays the role of the cut-off and, 
remarkably, for a hnite K the total number of states is hnite [^|. The continuum limit is 
recovered as L and K are sent to inhnity. For the sake of illustration, we display the set of 
states with K = 3. There are hve bosonic states 

V3^Tr[V(lM'(lMt(l)] |0> 

1 


and five fermionic states 


ii)» 

|2)i, 

|3)i, 

|4)» 

|5)i 

ID/ 

| 2 >/ 

|3>/ 

|4>/ 

| 8 >/ 


Ar3/2 

= —Tr 

N 

1 


Tr 


V(l)Bt(D-S'(l) |0> 


V(2)V(1)1 |0) 


-H' [b'(2)b'(D] |0> 

TV [rl'(3)] |0) 


iVl/2 


1 


Ar3/2 
1 


Tr [kl+(l)klt(l)pt(i)] |o) 

V37^Tr[Bt(l)Bt(l)Bt(l)] 10) 

[V(2)Bt(l)] |0> 


N 

= —Tr 

N 


'b\2)A\1)\ |0) 


AT/2 


Tr 


Pt(3)l |0). 


(4) 


For higher values of K all states may be generated systematically by considering all partitions 
of K into positive integers, assigning one of the two oscillators, or B\ to each integer, 
and eliminating null and redundant states. Once the basis of states for a given K has been 
found, it is straightforward to determine the matrix elements of I. 

As in the bosonic case, the mode operators create or annihilate “string bits,” except now 
there are two species of such bits. The stringy states we are interested in may be thought 
of as bound states of these string bits. Various terms in the supercharges implement joining 
and splitting of neighboring bits or an interchange of a bosonic and a fermionic bits. 

In what follows, we describe the result of our analysis for the spectrum of J = \/2iQ+Q-. 
We constructed the exact matrix representation of I symbolically for hnite K and then 
evaluated its spectrum numerically. The goal is to examine the convergence of the spectrum 
in the large K limit in order to recover the continuum theory. Unfortunately, the dimension 
of matrix I grows rapidly with K. Our computer resources have allowed us to perform 
diagonalizations for iF = 3, 4, 5, 6, 7, 8, 9, and 10, where the dimensions of the matrix I 
are 5, 10, 25, 62, 157, 410, 1097, and 2954, respectively. The number of states for each K 
is greater than that found in by one, since our model is not gauged and we include the 
single string bit state. 
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In figure we plot the spectrum of / as a function of y for K ranging form 5 to 10 
(we show only the eigenvalues whose magnitude does not exceed 3). Note that, for each 



.3 _^_I_ • - _I. . ' " _ ■ ..I _I_ iL— _ I . .. 

5 6 7 8 9 10 11 

Figure 1: Low lying spectrum of I for 0 < ?/ < 1 for K from 5 to 10 


K, there is a special value of y where the hrst zero eigenvalue appears. It is plausible 
that this phenomenon persists to arbitrarily large values of K, and that this zero crossing 
signihes a phase transition. Indeed, we found the locations of these zero crossings apparently 
converging like a power with increasing K. We have calculated these locations and tried a 
least squares £t to the functional form a + h/K^. This gives an estimate of yc = 0.37 for the 
zero crossing of the hrst level in the large K limit. 

We are also interested in the behavior of other states with low mass-squared. If a con¬ 
tinuum of states is forming as we conjectured, the n-th level for any hnite n should develop 
a zero eigenvalue aX y = yc in the limit K ^ 00 . This behavior is not at all obvious from 
hgure To examine the situation more closely, we found the eigenvalues of other low lying 
states for the value of y which corresponds to the zero crossing of lowest state and studied 
their dependence on K. These points are marked on hgure by tick marks on the spectrum. 
We demonstrate the K dependence of these quantities in hgure |[ Figure is plotted on 
a log-log scale in anticipation of the power law dependence. The linearity of the graph for 
the second level appears to confirm this. The third and fourth level also appear to follow 


7 

















Figure 2: Spectrum 2nd, 3rd, and 4th levels at the zero of 1st level v.s. K 


this trend superposed on some oscillations. The oscillation is likely due to the difference 
between the regularizations with odd and even K. For high enough K this oscillation should 
be suppressed and these low lying states will hopefully exhibit a power law decay similar to 
the second level. Alternatively we could extrapolate for even and odd K separately, but we 
do not have enough data for that. 

As a hnal check, we concentrated on the two lowest states and attempted to extrapolate 
their masses to inhnite K for each y. Again, the spectrum appeared to be converging like a 
power. Using best hts to the form a + 6/A''^, we determined our best guesses for the large K 
limit of the two lowest masses. The result of this calculation is illustrated in hgure where 
we plot the dependence on y of the two lowest masses (we show the results for hxed K as well 
as the extrapolations to inhnite K). We are encouraged to hnd that the extrapolated masses 
for the two lowest states come very close to vanishing simultaneously at y near 0.37, which 
is also the value of yc that we found by extrapolating the locations of the zero crossings. 

At the moment our results are far from being thoroughly convincing, since our extrapo¬ 
lations are no more than the best guesses we can offer based on the available computational 
data. With the computer resources available to us we were unable to extend the calculation 
to K > 10. Bigger computers, or improved algorithms, should certainly permit the extension 





Figure 3: Spectrum of two lightest states for various iF’s and their extrapolation. Dotted 
lines are for K from 3 to 10. The solid line is the extrapolation. 


of our calculations to higher K. It would certainly be interesting to check the continuation 
of the trends we observed and perform more reliable extrapolations. 

In conclusion, we have presented some numerical evidence which suggests that the spec¬ 
trum of a supersymmetric large-iV matrix held theory becomes continuous and massless at 
a critical value of the coupling constant, Ac. A similar phenomenon occurs in the c = 1 
matrix model where Ac corresponds to the onset of the world sheet continuum limit, and the 
Liouvihe mode emerges as an additional target space coordinate. By analogy, we speculate 
that our H-l dimensional matrix model describes a 2-1-1 dimensional non-critical superstring 
in some background. If true, one might further consider an exciting possibility of hnding a 
non-perturbative matrix model formulation of higher dimensional superstring theories. 

By construction, the conjectured 2-|-l theory will possess (1,1) supersymmetry in two of 
its dimensions. From the special nature of the Liouvihe mode, however, we hnd it extremely 
unlikely that this symmetry is extended to 2-1-1 dimensions. As such, it is similar in spirit 
to the model considered in |^. It would be interesting to further explore the relationship 
between these two models. Perhaps more can be learned about the relation between our 
approach and its world sheet formulation by studying the renormalization group flow and 
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the operator algebra of the world sheet action (|^ along the lines of []^. It would also 
be interesting to explore the relation to the Green-Schwartz approach taken in P, The 
NSR approach, the Green-Schwartz approach, and the matrix models offer complementary 
methods for the dehnition of supersymmetric non-critical strings. It would be very interesting 
to hnd a model which we understand in all three formulations. 
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